Super-Eddington accreting massive black holes as long-lived cosmological standards 
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Super-Eddington accreting massive black holes (SEAMBHs) reach saturated luminosities above a certain 
accretion rate due to photon trapping and advection in slim accretion disks. We show that these SEAMBHs 
could provide a new tool for estimating cosmological distances if they are properly identified by hard X-ray 
observations, in particular by the slope of their 2-10 keV continuum. To verify this idea we obtained black 
hole mass estimates and X-ray data for a sample of 60 narrow line Seyfert 1 galaxies that we consider to be the 
most promising SEAMBH candidates. We demonstrate that the distances derived by the new method for the 
objects in the sample get closer to the standard luminosity distances as the hard X-ray continuum gets steeper. 
The results allow us to analyze the requirements for using the method in future samples of active black holes 
and to demonstrate that the expected uncertainty, given large enough samples, can make them into a useful, new 
cosmological ruler. 
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The discovery of the accelerating expansion of the Universe 
has now been established through observations of type la su- 
pernovae (SNe la) 01, and is likely to be confirmed further 
with new standard rulers provided by baryon acoustic oscil- 
lations, weak lensing and clusters of galaxies J2]. However, 
SNe la beyond z > 1.5 are rare as there is no time for 
their progenitors to evolve in substantial numbers given their 
lower metallicity J3l ■ To further probe the dynamics of the ac- 
celeration, new distance indicators are needed at and beyond 
these redshifts. Based on well-understood physics we show in 
this Letter that super-Eddington accreting massive black holes 
(hereafter SEAMBHs) in some active galactic nuclei (AGNs), 
that are characterized by a mass of 10 6 ~ 8 Mq, can provide a 
new tool to estimate cosmological distances at a wide range 
of redshifts, including the high redshift Universe. 

Radiation pressure limits the spherical accretion rate onto 
black holes to MBdd = -^Edd/^c 2 , where Leaa = 
4irGM,m p c/cr T is the Eddington luminosity for a pure hy- 
drogen plasma, rj (~ 0.1) is the mass to radiation conversion 
efficiency, er T is the Thomson cross section, m p is the proton 
mass, c is the speed of light, G is the gravitational constant 
and M, is the black hole mass. However, super-Eddington 
accretion onto black holes is feasible in slim disks where the 
radiation pressure-dominated regions (RPDR) are thermally 
stable due to the radial advection of the locally emitted ra- 
diation Qj. In such disks, the timescale of photon diffusion 
to the disk surface is longer than that of the radial motion of 
the accreting gas in the RPDR. The photons are trapped in- 
side the accretion flows and are advected into the black holes. 
This advection dominates within the photon trapping radius, 
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rather than linearly) to the accretion rate [0 0], 
L. = £ (l + olnmi 5 )M. , 



(1) 



Af. is the mass accretion rate and R g = GA/./c (6fl. Photon 
trapping affects the total emitted radiation and results in a sat- 
urated luminosity, L,, which is proportional (logarithmically, 



whereto « 5.29 xlO 38 ergs" 1 Mq 1 , and ops 0.476 0]. For 
reference, at m = 15 the saturated luminosity is ~ 4.20LEdd- 
Thus, at a given black hole mass, SEAMBHs are radiating 
basically at a constant luminosity which, as shown below, can 
therefore be used to deduce cosmological distances. 

In this Letter we address two important issues: how to iden- 
tify SEAMBHs and how to test, observationally, Equation [TJ 
and its uncertainties such that it can be used to derive reli- 
able cosmological distances. While SEAMBHs arepredicted 
to have unique optical-UV spectral characteristics [ , their 
use to identify such sources is hampered by the dilution of 
the disk emission by stellar radiation from the host galaxy at 
long wavelengths, and by the Galactic and inter-galactic ab- 
sorption at short wavelengths. In fact, current observations 
cannot identify such systems using only their spectra in the 
optical-UV domain. Fortunately, X-ray spectroscopy allows 
such identifications for two reasons. First, there is a well- 
known positive correlation between the 2-10 keV X-ray pho- 
ton spectral index (T) and the Eddington ratio (Lboi/ -^Edd) 
IH]. In addition, higher LBoi/^Edd sources emit a smaller 
fraction of their total radiation (/ x = Lx/L,) at hard X-ray 
energies These properties are easy to measure with mod- 
ern X-ray observations and are similar to those observed in 
Galactic black holes in X-ray binaries ifioll . 

The general theory that links X-ray emission to the optical- 
UV spectrum of accretion disks is based on the assumption of 
a hot corona above the disk that is the source of the X-ray ra- 
diation. The X-ray emission efficiency of the processes in the 
corona depends on ieoi/^Edd- ^ n particular, the magneto- 
rotational instability is a key factor to produce the viscosity 
that helps transporting angular momentum outward lUll . The 
process takes place through magnetic buoyant transportation 
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above the cold disk 11211 which leads to a "corona-dominated" 
dissipation through hard X-ray emission. The increases of the 
accretion rates result in the weakening of the transportation of 
magnetic tubes because of the inflation of the disk by radia- 
tion pressure which yields a reduction of the buoyant velocity. 
The end result of these processes is the radial advection of the 
emitted photons and the suppression of the relative X-ray flux 
(fx) °f the system 111 311 . 

The second effect of an increasing Leoi/^Edd ratio is the 
steepening of the X-ray photon index V. There have been var- 
ious attempts at calculating V in hot coronae from basic con- 
siderations of the conditions in accretion disks 111 211 . The 2-10 
keV emission is mainly due to the Comptonization of photons 
from the cold disks (the X-ray reflection can be neglected in 



this band 1 14]). It has been shown that / x oc M, 1 ^ 18 m -4 / 9 
in the RPDR (see equation 13 for / x < 1 and Figure 1 in Ref. 
1 1311) and thus one expects that a luminosity f K L, is radiated 
by the coronae. Since Comptonization is the main cooling 
process, the balance between heating and cooling yields the 
density of hot electrons n c oc / x and n c oc m~ 4 / 9 . 

Under the conditions discussed above, the X-ray photon 
index can be approximated by V ~ 2.25y _2/ ' 9 IU5I1 . where 
y = 4# e r T (l + 4# e )(l + t t ) is the Comptonization param- 
eter, 8 e = kT e /m e c 2 , T e is the electron temperature and r T 
is the Thomson scattering optical depth. It is expected that 
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where 7 = 1 for unsaturated Comp- 
tonization (r T < 1) and 7 = 2 for saturated Comptonization. 
This scenario is supported by the behaviors of black hole X- 
ray binaries in very high states [[loll , showing that the 2-10 
keV emission is at the level of low/hard states, but the pho- 
ton indexes are typically T > 2. Considering that hard X- 
ray spectra have a cutoff of 50-100 keV (0 e - 0.1 - 0.2) 
IU6I1 and applying the above coronal model to AGNs, we 
have r oc m 87 / 81 . This prediction agrees with the observed 
r — LBoi/^Edd correlation of AGNs hosting standard (opti- 
cally thick and geometrically thin) accretion disks ]8|]. 

The simplest model for the slim disks of SEAMBHs 
assumes a spherical hot corona with a characteristic size 
£ c . The Comptonization of photons from the slim 
disk surface produces a hard X-ray luminosity Lx = 
A6 e n e cr T c (Z./4x^c) (47r^/3), giving rise to t t = 
3/x/46» e w 0.8 fo.i&o.l, where r T = n e a T £ c , / .i = / x /0.1 
and #0.1 = e /0.1. The coronae having a Comptonization 
parameter of y ps 0.8, the SEAMBHs are then characterized 
by T > 2.3y g^ 9 , where y„ 8 = y/0.8. Obviously there are 
uncertainties in these parameters and the resulting theoretical 
relationships, however, we can use the observed T to identify 
SEAMBHs in different types of active galactic nuclei. 

The best group of AGNs where such processes have been 
studied are narrow line Seyfert 1 galaxies (NLSls). These ob- 
jects are separated from broad line Seyfert 1 galaxies (BLSls) 
by having the following properties lil7[ll8ll : 1) the full-width- 
half-maximum (FWHM) of H/3 profiles are narrower than 
2000 km s _1 ; 2) strong soft X-ray excess; 3) unusually strong 
(relative to H/3) optical iron emission lines; 4) weak [O III] 



lines ([Olll]/H/3 < 3); 5) fast, large amplitude X-ray vari- 
ations lfl9ll . The typical black hole mass in NLSls (see be- 
low) is considerably smaller than that in BLSls of similar 
Lbo\ which implies that many of them may be accreting at 
super- Eddington rates. 

To test observationally EquationQ]in NLSls, we have to es- 
timate M,. Black hole masses can be measured individually 
using the reverberation mapping (RM) technique, which in- 
vokes the response (time lag) of the broad emission lines with 
respect to changes in the continuum produced by the under- 
lying disk |2(| 21]. In its most detailed version, the velocity- 
resolved reverberation mapping (VRRM), one can derive the 
spatial distribution of the line emitting gas, and its velocity, at 
every location around the black hole. In principle, this phase- 
space mapping enables us to determine, accurately, the black 
hole mass (this is equivalent to dynamical methods which are 
used to measure the masses of black holes in normal nearby 
galaxies). Such two dimensional mappings are only available 
for a handful of sub-Eddington sources and many details of 
the technique need to be improved for accurate measurements 
(e.g. the case of the BLS1 Mrk 50 JH). 

Fortunately, we can use as a proxy the tight correlation be- 
tween the size of the broad line regions (BLR) (the time lag 
times the speed of light) and the underlying continuum lu- 
minosity to obtain an empirical relationship, for a sample of 
about 35 AGNs, that can be combined with the observed line 
widths to estimate the black hole masses in a large sample of 
sources. The relationship is given by 
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10 44 erg s" 1 



(2) 



where R ~ 9 x 10 °cm, a = 0.6 ± 0.1, and L5100 is 
the AGN continuum luminosity (XL\ at 5100 A in units of 
10 44 erg s" 1 ) H IH. We note that among the 35 AGNs 



used to derive this correlations, eight are NLSls and follow 
the same trend as other sources |]23J. 

Equation|2]enables us to obtain the black hole masses by as- 
suming a virialized (gravitationally bound) cloud system and 
combining all the unknown geometrical factors, such as the 
inclination to the line of sight, into a single constant / BLR . 
Using this constant we can now write an expression for the 
"reverberation mapping-based virial mass estimate" of super- 
massive black holes as 



M. = /b^G" 1 ^ 



(3) 



where U FW hm is the FWHM of the broad emission line (e.g. 
H/3) that was used to derive the time lag in the RM measure- 
ment. The factor / BLR is calibrated by comparing the results of 
the RM experiments to direct black hole measurements based 
on the M. — 0* relation, where er» is the stellar velocity dis- 
persion in the bulge of the host galaxy. Such a comparison 
is now available for about 30 out of the 35 AGNs in the RM 
sample. It shows that / BLR ~ 1.2 ± 0.2 J3. Noting that 
L, = k b L^ioo, where k b is a bolometric correction factor, 
we obtain an expression of L, in terms of / BLR and Vf WH m- 
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FIG. 1: Left: Residuals of the SEAMBH distance modulus (Afi = fi, — /iL,). The error bai's are taken as A/i, = 1.17 mag (see text for details). The assumed 
cosmological parameters are Ho = 71 km s" 1 Mpc -1 , Ca = 0.73 and !)m = 0.27. Middle: The distribution of Afi as a function of the hard X-ray (2-10 
keV) photon index I\ The normalized frequency is referred to fractions to the peak number of objects. N is the total number of SEAMBHs selected by I\ 
Right: scatter of the Afi distributions with V (bottom) and the distribution of T for the sample of SEAMBH candidates. The dependence of the dispersion in the 
residuals as a function of V shows a systematic decrease and tends to cr^™, indicating the efficacy of the selection as standard candles. 



Since L5100 = 47r2? 2 i 7 5ioo, where -F5100 is the measured con- 
tinuum flux in units of erg s" 1 cm -2 at 5100 A and T> % is the 
luminosity distance of the black hole, we obtain the expres- 
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This expression still involves the unknown accretion rate rh.15 
which, as suggested earlier, can be estimated from the X-ray 
slope r. However, the dependence for T > 2 (which is the one 
we are interested here, see below) is weak enough that we can 
use the approximation TO15 = 1 and obtain P. from the other 
known observables and constants. Thus, having established 
that an object is a SEAMBH by measuring its T index, we 
have a way to directly measure its distance. 

In the following we use the distance modulus, fi. = 
51og(2?./pc) — 5, and compare it with the one obtained 
from the standard luminosity distance 2?l in the Friedman- 
Lemaitre-Robertson- Walker metric /i^ = 51og(2?L/pc) — 5. 
The prediction is that comparing V, and L>l we will get 
smaller residuals Afi = fi, — fi^ f° r larger T, since large 
indices point to conditions closer to those predicted by the 
SEAMBH theory. 

We now turn to the available samples of SEAMBH can- 
didates. While the observed NLS1 properties may all be re- 
lated to the large Eddington ratio 11181 12511 . not all NLSls have 
super-Eddington accretion rates. It is thus necessary to use 
the hard X-ray spectra to identify SEAMBHs among NLSls. 
We selected a lar ge n umber of NLSls from several hetero- 
geneous samples |18, 26| with hard X-ray observations by 
ASCA, XMM-Newton, Chandra and Swift \T&. All data and 
data reduction details are provided in the Supplementary Ma- 
terials Online. In short, we use the observed .F5100 flux and 
the estimated black hole mass M, to derive n B and hence L, 
for each source (see also a somewhat different approach but 



similar results in 17|,|28|]). The contaminations by the stellar 
light of the host galaxy were removed, prior to the estimate 
of -F5100, using the approximation described in Materials On- 
line. We examined the LBoi/^Edd distribution in our NLS1 
sample and found that, indeed, many of them indicate super- 
Eddington ratios, up to 5 and even more. Thus, the selection 
of sources by T is indeed a good way to identify such sources. 

We calculate Afi for all sources, bin them into various 
groups of different T, and plot them in FigureQ] As shown by 
the standard deviation (ca^X me scatter of Afi systematically 
decreases as V increases. From Equation|4]i, this behavior can 
be understood in terms of the scatter expressed by 
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where A(i x = J^Li = A (AXJX. 

/blr, Rq, t B , V FvmM , ^5100. A(i = 1,2,3) = 1/2(1 - a), 
A4 = 1/(1 — a) and A$ = 1/2. A^. converges to Afix as 
Am/m decreases with increasing T. 

To illustrate the typical uncertainty on individual points we 
assume that A/ BLR // BLR = 0.2 12411 (from the scatter in the 
M, — a* distribution), AF5100 /-F5100 = 0.2 (from the known 
variations in the optical continuum and the uncertainty in the 
substraction of the stellar background), AT4- WHM /Vi. W HM = 
0.05 (from the uncertainty in fitting the emission line profiles 
and measuring V FWHM ), AR /R = 0.2 and Ak b /k b = 0.3. 
The estimated uncertainty of k b is the largest and most prob- 
lematic for several reasons. First, all our estimates of k b are 
derived from theoretical calculations of slim disk spectra 10]. 
These have not been verified observationally because of the 
lack of extreme UV observations, where most of the emitted 
luminosity in disks around small black holes (10 6 — 10 7 M Q ) 
is emitted. Moreover, a factor of 10 increase in black hole 
mass results in a factor of IO 1 / 3 decrease in k b (from about 
100 in 10 6 M Q to about 40 in 10 7 M Q black holes). The entire 
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range of black hole masses in our sample suggests a very large 
Ak b /k b . Fortunately, the individual masses are known and 
the uncertainty on massive black holes using the RM-based 
virial method is only a factor of about 3 112911 - This and the 
allowed range of spectral shapes and tin gives the quoted esti- 
mate on Ak b /k b . 

Combining all these uncertainties and assuming a = 0.6 
we obtain Afix = 0.54. This corresponds to Afi, ~ 1.17 
mag since the first term in Equation [TJ tends to vanish for 
large enough values of T. This is not surprising given Equa- 
tion Q] and the known uncertainties on AI, measured by the 
RM-based virial method. This uncertainty is estimated to be 
0.3-0.5 dex which would suggest a similar uncertainty on L m . 
Since fi oc 2.5 log L,, we get a similar uncertainty on individ- 
ually measured points to the one obtained by the more detailed 
calculations. From Figure [TJ we find A/i. > a^™ = 0.93 
mag, which clearly shows that the selected sample has a small 
scatter of Arn/m related to this term. 

The panels of Figure [TJ illustrate both the convergence to 
(A/i) = and the reduced scatter when using increasinglt 
larger values of V. The central panel shows that for a suffi- 
ciently steep X-ray continuum, a combination of a large num- 
ber of SEAMBHs gives, indeed, the correct distance with a 
small scatter. 

The systematic decrease of the dispersion in the subsam- 
ples with increasing T, while keeping a median with little 
variation, cannot be accounted by statistical fluctuations in 
the subsamples. Considering that the subsamples are not in- 
dependent, we assess the statistical significance using Monte 
Carlo simulations with 10 7 Gaussian samples of size N = 60 
with the same mean and dispersion as the observed one. From 
them, we selected distinct (i.e., with no replacement) random 
sub-samples of the same size as the ones selected (Figure [TJ, 
and estimate the probabilities that the differences in the medi- 
ans and the ratio of the dispersions correspond to the observed 
ones. We find that while the probability that the medians do 
not differ from the one of the underlying sample of N = 60 
is always very large (above ~ 80%), the probability that the 
dispersions are as small as the observed one is always smaller 
than 3%. The observed trend cannot therefore be ascribed to 
random fluctuations of small samples extracted from the main 
sample. In the sample used here, there are only 12 sources 
with r e [2.3, 2.5] but future samples will be larger since 
such objects can be observed to high redshift. We point out 
that the T > 2.3 and 2.3 < T < 2.5 panels are statistically 
indistinguishable because of the poor quality of the sources 
with T > 2.5. There are 11 sources in total with T > 2.5 
listed in Table 1 in the Supplementary Material Online. Five 
of them with hard X-ray observations have large error bars 
(Ar > 0.37), making the T— binning less significant for these 
small samples. Five other objects are observed in the 0.5-8 
keV band by Chandra (Williams et al. 2004 in II2611 . see notes 
in Table 1), but the data quality only allows us to approximate 
the 0.5-8 keV spectral indexes with the 2-8 keV proxy, and 
the error bars of these sources remain very large (AT > 0.3 
except one). The 2.3 < T < 2.5 panel is shown for refer- 



ence and future observations. Moreover, the accuracy of the 
measured M, can increase substantially if new, dedicated RM 
experiments are carried out on a large number of SEAMBH 
candidates. This can reduce the uncertainty on / BLR , -F5100, 
Vp WHM and Rq and constrain a to be the slope for this popu- 
lation only. As FigureQJshows, it is reasonable to assume that 
for a large enough sample of SEAMBHs at a narrow redshift 
range, we could expect a scatter on the population /.t. that ap- 
proaches 0.15 mag, similar to the current accuracy of SNe la 
method Jl. 

SEAMBHs, as a new type of cosmological distance indica- 
tor, have a number of advantages over others 113 Oil : (1) The sat- 
urated luminosities are well understood on physical grounds 
and have no potential cosmic evolution. (2) The objects can 
be efficiently selected by their X-ray or optical properties. (3) 
They can probe a large range of redshifts, as they follow the 
cosmic growth of massive black holes that are abundant at 
high-z and are very luminous. (4) Unlike SNe la, repeated 
observations can be made to improve the observational accu- 
racy. 

Several observational issues require careful attention: (1) 
As mentioned before, the R BLR — L relation applied here is 
the one obtained for all AGNs. The re-calibration of this re- 
lationship in a dedicated NLSls or SEAMBH sample, by ob- 
taining better estimates of Rq and a in Equation|2] can reduce 
the scatter in the derived mass and hence k b . Obviously, the 
estimated M. involves the distance to the source which is the 
quantity we are attempting to measure. However, we do not 
expect large differences between 2?l and T>, so this uncer- 
tainty, by itself, is very small. (2) The 2-10 keV luminosity 
and slope are both variable IU9I1 which may lead to misiden- 
tification of SEAMBHs. Long-term averaged values can be 
used to improve the accuracy. 

The prospects of building large samples of SEAMBHs to 
be used as tests of the current cosmological model are promis- 
ing. We expect that roughly 20% of NLS1 with T > 2 host 
SEAMBHs. As NLS1 constitute about 10% of all AGNs, 
there should be 10 4 ~ 10 5 SEAMBHs among Seyfert 1 galax- 
ies with z < 0.3 13 ill . SEAMBHs could be even more abun- 



dant at high-z although the definition of NLS 1 should be mod- 
ified in such cases 113211 . Here we require black hole mass esti- 
mates that are based on both the H/3 and Mg II A2798 A lines. 
X-ray spectra can be obtained by Nustar, by the upcoming 
eRosita and HXMT missions 113311 . Given accurate observa- 
tions of SEAMBHs at high redshift we will have a unique 
chance to explore in-depth the dynamics of the accelerating 
Universe as well as the nature of dark energy. 
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